The effects of nuclear spins on the quantum relaxation of the magnetization 

for the molecular nanomagnet Fcg 
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The strong influence of nuclear spins on resonant quantum tunneling in the molecular cluster Feg 
is demonstrated for the first time by comparing the relaxation rate of the standard Feg sample with 
two isotopic modified samples: (i) ^®Fe is replaced by ^'^Fe, and (ii) a fraction of is replaced 
by ^H. By using a recently developed "hole digging" method, we measured an intrinsic broadening 
which is driven by the hyperfine fields. Our measurements are in good agreement with numerical 
hyperfine calculations. For T > 1.5 K, the influence of nuclear spins on the relaxation rate is less 
important, suggesting that spin-phonon coupling dominates the relaxation rate. 



PACS numbers: 75.45.+j, 75.60Ej 



Mesoscopic quantum phenomena are actively investi- 
gated both for fundamental science and for future appli- 
cations, for instance in quantum computing. Magnetic 
molecular clusters are among the most promising candi- 
dates to observe mesoscopic quantum phenomena 
One of the most prominent examples is an octanuclear 
iron(III) cluster, called Fes (Fig. [^), with a spin ground 
state of S" = 10 Q|. Below 360 mK, the magnetization 
relaxes through a pure tunneling process giving rise to 
a stepped hysteresis cycle [Q. Furthermore, the tun- 
nel splitting A of Fes shows periodic oscillations when 
a transverse magnetic field is applied along the hard 
axis a long searched phenomenon in magnetism as- 
sociated with the Berry phase ||], and predicted sev- 
eral years before [Q. Since A is extremely small for the 
ground state tunneling, ca. 10^^ K at = 0, the tun- 
neling process should occur only in an extremely narrow 
magnetic field range, ca. 10^^ T, and should be prac- 
tically unobservable. However, a recent theory proposes 
that the tunneling is mediated by fluctuating hyperfine 
fields generated by magnetic nuclei but direct exper- 
imental evidence is so far lacking. 

In order to study the influence of nuclear spins, we in- 
creased the hyperfine coupling by the substitution of ^^Fe 
with ^^Fe, and decreased it by the substitution of with 
^H. We found that the relaxation rate of magnetization 
in the tunneling regime shows a clear isotope effect which 
we attribute to the changed hyperfine coupling. 

The crystals of the standard Fe8 cluster, '^'Fes or Fes, 
[Fe8(tacn)602(OH)i2]Br8.9H20 where tacn = 1,4,7- tri- 
azacyclononane, were prepared as reported by Wieghardt 
et al. . For the synthesis of the ^^Fe-enriched sample, 
^^Fes, a 13 mg foil of 95% inriched ^^Fe was dissolved in a 
few drops of HCI/HNO3 (3:1) and the resulting solution 
was used as the iron source in the standard procedure. 



The ^H-enriched Fes sample, °Fes, was crystallized from 
pyridine-ds and D2O (99%) under an inert atmosphere 
at 5°C by using a non-deuterated Fe(tacn)Cl3 precursor. 
The amount of isotope exchange was not quantitatively 
evaluated, but it can be reasonably assumed that the H 
atoms of H2O and of the bridging OH groups, as well as 
a part of those of the NH groups of the tacn ligands are 
replaced by deuterium while the aliphatic hydrogens are 
essentially not affected. The crystalline materials were 
carefully checked by elemental analysis and single-crystal 
X-ray diffraction. 




FIG. 1. Schematic view of the magnetic core of the Feg 
cluster. The oxygen atoms are black, the nitrogen gray and 
carbon atoms are white. For the sake of clarity only the hy- 
drogen atoms that are exchanged with deuterium are shown as 
small spheres. The arrows represent the spin structure of the 
ground state 5* = 10 as experimentally determined through 
polarized neutron diffraction experiments . 

The magnetic measurements were made on single- 
crystal samples by using an array of micro-SQUIDs [ p^ , 
which measure the magnetic field induced by the magne- 
tization of the crystal. The advantage of this magnctome- 
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ter lies mainly in its high sensitivity and fast response, 
allowing short-time measurements down to 1 ms. Fur- 
thermore the magnetic field can be changed rapidly and 
along any direction. 
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FIG. 2. Comparison of the relaxation rates of three differ- 
ent Feg samples. The time needed to relax one percent of 
the saturation magnetisation Ms is plotted versus inverse tem- 
perature 1/T. The initial magnetisation Minit was reached by 
a fast cooling in zero applied field. The relaxation was mea- 
sured in a field of ^oH^ — 42 mT. In the inset typical relax- 
ation curves (M versus time) recorded for the three sample 
at T = 40 mK are shown. 
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FIG. 3. Comparison of the short time relaxation rates of 
three different Fes samples at T = 40 mK with Htrans ~ and 
Minit ~ 0. The inset displays an typical example of a hole 
which was dug into the distribution by allowing the sample 
to relax for the time t^ig at fioH^ig ~ 14 mT. 

In order to avoid the influence of the crystal shape ||Tl| ] 
the relaxation was measured starting at an initial mag- 
netization Afinit = where intermolecular dipolar in- 
teractions lead to a field distribution with a width of 
about 50 mT [|lO| . A small field Hz was then applied and 
the relaxation of magnetization was measured (inset of 
Fig. H). For all three samples the relaxation was clearly 
non-exponential and could be adjusted to a -^/Fsqitt— law 
for the short time regime (t < 100 s) for T < 0.4 K. Fig. ^ 
displays the field dependance of T sqit (Hz). The relax- 



ation of the three samples at 40 mK are strikingly differ- 
ent from each other. The ^^Fes sample is the fastest re- 
laxing one whereas the °Feg shows the slowest relaxation 
rate. As a complete theory of the relaxation behavior of 
crystals of molecular clusters is still missing ||l2| , we plot 
in Fig. Ijthe time needed to relax one percent of the satu- 
ration magnetization Mg as a function of the inverse tem- 
perature 1/T. Relaxation and ac susceptibility measure- 
ments at r > 1.5 K showed no clear difference between 
the three samples suggesting that above this tempera- 
ture the relaxation is predominately due to spin-phonon 
coupling ||l^,|ljl . Although the increased mass of the iso- 
topes changes the spin-phonon coupling, we believe that 
this effect is small. 

In principle, the change of masse does not change the 
crystalline field of the Fe ions, i.e. the anisotropy con- 
stants. Experimentally, this is confirmed with measure- 
ments below T < 0.35 K, where spin-phonon coupling 
is negligible, by two observations: (i) relative positions 
of the resonances as a function of the longitudinal field 
Hz are unchanged ||l^ , and (ii) all three samples showed 
the same period of oscillation of A as a function of the 
transverse field Hx , a period which is very sensitive to 
any change of the anisotropy constants. Finally, we point 
out that the mass is increased in both isotopically modi- 
fied samples whereas the effect on the relaxation rate is 
opposite. 

A deeper insight into the relaxation mechanism can 
be achieved by using our recently developed hole digging 
method which allows us to estimate the hyperfine level 
broadening ||l^. Starting from the well defined magne- 
tization state Minit = 0, and after applying a small field 
Hdig, the sample is let to relax for a time t^ig, called 
digging field and digging time, respectively. During the 
digging time, a small fraction AMdig of the molecular 
spins tunnel and reverse the direction of their magneti- 
zation. Finally, a field H is applied to measure the short 
time square root relaxation rate Fsqi-t [ p^|JTl|] . The entire 
procedure is then repeated to probe the distribution at 
other fields yielding the field dependence of the relaxation 
rate TsqrtiH, iJdig, ^dig) which is more or less proportional 
to the number of spins which are still free for tunneling. 
The result of this procedure is that a very sharp "hole" 
is dug into the rather broad distribution of Fsqrt |lO) . A 
typical example is shown in the inset of Fig. ^. 

In the limit of very short digging times, the differ- 
ence between the relaxation rate in the absence and in 
the presence of digging, Fhoic = T^sqrtiH, Hdig,tdig = 
0) — Fsqi.t(i?, -ffdig, ^dig), is approximately proportional to 
the number of molecules which reversed their magneti- 
zation during the time tjig. Fhoie is characterized by a 
width that wc call the hole line width a. In order to 
find a hole line width that is close to the hyperfine level 
broadening, all the effects that can broaden the measured 
hole width must be reduced. The experimental condition 
giving the smallest line width was found for hole dig- 
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ging in the tails of the dipolar distribution and for small 
initial magnetization. Under these conditions the spins 
that tunnel are statistically far from each other allowing 
us to measure tunneling in the diluted limit. In addi- 
tion, we applied a transverse field of /J.o-fftrans = 200 niT 
parallel to the hard axis which reduces the tunnel rate 
allowing us to dig very tiny holes. Fig. ^ displays the 
hole line width cr as a function of the reversed fraction 
AMdig/2Ms of molecular spins. A linear extrapolation of 
(7 to AMdig/2Ms = gives cto which is directly associated 
to the hyperfine level broadening |^ . Experimentally we 
found ao to be 0.6 ± 0.1, 0.8 ± 0.1, and 1.2 ± 0.1 mT, for 
^Fes, ''^Fes, and ^''Fes, respectively. 
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FIG. 4. Hole line width cr as a function of reversed fraction 
AMdig to dig the hole. The initial thermal distribution of 
diploar fields was reached by a fast cooling in zero applied 
field. In order to reduce the line width the hole was dug at 
a longitudinal field of /io-ffz = 42 mT, and in the presence 
of a transverse field of /lo-fftrans = 200 mT applied along the 
hard axis. In the inset shows typical hole shapes observed at 
AMdig = 0.001 Ms 

The isotope effect, observed here for the first time in 
magnetic nanostructures, clearly point out the role of the 
magnetic nuclei in the relaxation of the magnetization. 
An evaluation of the hyperfine fields in the three different 
samples is therefore necessary. The hyperfine interaction 
between the total spin S of the cluster and the magnetic 
nuclei can be decomposed into the sum of terms related 
to the interactions between the magnetic moment li of 
the I— th nucleus and the individual spin 5^-, assumed to 
be localized on the j— th iron centre flq]: 



Bn.f - = S E E ^^^^^ 



(1) 



The projection coefficients Cj in Eq. (|I|) depend on the 
wave function of the ground state S = 10, which can 
be calculated by diagonalizing the S* = 10 block (6328 x 
6328) of the exchange spin-Hamiltonian matrix oi H = 



With the exchange-coupling parameters that best re- 
produce the temperature dependence of the magnetic 
susceptibility |]l^, the spin configuration depicted in 
Fig. 1^ provides a large contribution (of ca. 70%) of 
the ground state S = 10. According to this picture, 
which has been recently confirmed by polarized neu- 
tron diffraction data ||l^, the projection coefficients are 
C3 = C4 = —5/22 for the spins pointing down and 
ci = C2 = = ... — cs — 8/33 for the remaining iron 
spins. 

The hyperfine interaction described by the Ay tensors 
is both through-space (dipolar) and through bond (con- 
tact) in nature. The dipolar components can be easily 
calculated with the point dipolar approximation. The 
most important coupling with the nuclei are those of 
the bridging OH groups having A^ constants as large as 
0.045 mT, while the hyperfine coupling with N and Br 
nuclei does not exceed 0.005 and 0.003 mT respectively. 
An order of magnitude for the contact terms was esti- 
mated using a Density- Functional Theory calculation |lj] 
at the B3LYP level ||l^ on a model symmetric dimer 
[(NH3)4Fe(OH)2Fe(NH3)4]^+ • 

The hyperfine interaction generates a field which splits 
the ms = ±10 states. Since each nuclear spin I splits 
each state into 21+1 sublevels, the number of sublevels 
generated by the coupling of 18 ^"'N atoms (/ = 1), 8 
™'8iBr atoms (/ = 3/2), and 120 atoms (I ==1/2) 
present in ''^Feg is prohibitively large, being 3^^x4^x2^^''. 
We therefore made an approximation taking into account 
only the most significant terms. In fact the total line 
width goes as the geometric sum of the individual contri- 
butions and therefore the largest contributions dominate 
over the smaller ones. We evaluated the gaussian broad- 
ening determined by the 12 nuclei of the OH groups, 
assuming them equivalent with a contact hyperfine cou- 
pling constant ^cont — 0.05 mT, of the 18 ^H nuclei of 
the NH groups assuming Acont — 0.025 mT, and of the 14 
N nuclei with Acont = 0.2 mT. These values introduced 
in Eq. (|]) gave gaussian lines with widths at half-height 
of 0.2, 0.15, and 0.4 mT, respectively. By combining 
these we estimate a resulting gaussian distribution with 
a line width of 0.5 mT, in acceptable agreement with the 
experimental value of 0.8 mT. 

The effect of the ^^Fe nuclear spins in the enriched 
samples was estimated by assuming that each nucleus 
only feels its own electron spin. Therefore Eq. (^) is 
simplified as only the terms with i — j are different from 
zero. Using ^('"'^Fe) = 1.0 mT, in agreement with re- 
ported data 1 20 1 , we calculate the stick diagram reported 
in Fig. ^ which arises from the coupling of six equiv- 
alent 5^Fe / = 1/2 with A = 8/33x1.0 mT and two 
other nuclear spins 7=1/2 with coupling A — — 5/22x 
1.0 mT. The related histogram can be fitted with a gaus- 
sian whose line width is ca 0.8 mT. If we consider the 
experimental line width of the resonance, ca. 0.8 mT, 
of '^'Fes and we add the contribution of the ^^Fe nuclei 
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we obtain AiJ « 1.1 mT in close agreement with the 
experimentally observed value of 1.2 mT. 

The partial substitution for the nuclei of the OH 
and NH groups with the less magnetic isotopes leads 
to a reduction of the line width which in our calculations 
is estimated to be ca. 0.1 mT which should be compared 
with the experimental narrowing of ca. 0.25 mT. The dif- 
ference between the calculated and observed reduction of 
the line width is similar to the smaller calculated line 
width compared to the observed one of the '^'Feg sample. 
This may come from an underestimation of the interac- 
tions with the nuclei. 

The present data show the fundamental role of the nu- 
clear spins in the relaxation of the magnetization of Fcg 
in the quantum regime. Indeed, this is in contrast to the 
familiar role of isotope substitution which are generally 
associated with phonon coupling and thus proportional 
to the mass of the nuclei. Here we show that it is the 
magnetic moment of the nuclei which is important at 
temperatures well above those at which nuclear spin po- 
larisation is observed. 

We are indebted to A. Bencini for his help in DFT 
calculations, and to P.C.E. Stamp, I. Tupitsyn, N.V. 
Prokof'ev, and J. Villain for many fruitful and motivat- 
ing discussions. The financial contributions of Italian 
MURST, CNR, PFMSTAII, and the French DRET are 
acknowledged. 
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FIG. 5. Calculated histogram of level splitting due to the 
coupling with ^^Fe 1 — 1/2 nuclear spins in the Feg clusters. 
In equation (jl|) Aij constant have been assumed to be 1.0 mT 
for i — j and zero for i ^ j. The solid line represents the best 
fit using a gaussian line 
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